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The linear response theory is used to describe magnetoresistance oscillations of short-period 
unilateral superlattices with strong modulation (or alternatively arrays of coupled quantum wires). 
The semiclassical description of this system fails for strong magnetic fields (magnetic breakdown) 
and we employ a simple fully-quantum-mechanical tight-binding model in conjunction with Kubo's 
formula instead. The resulting magnetoresistance data nicely compare to the experiments while the 
model opens good intuitive insight into the effects taking place in the system. 

PACS numbers: 73.63.Nm 



I. INTRODUCTION 

Transport properties of two-dimensional electron sys- 
tems (2DES) with unidirectional periodical modulation 
have been studied for more than ten years. Lot of at- 
tention has been paid to the case of weak modulation by 
electric and magnetic fields. Such a system was first pre= 
pared by means of holographic techniques by Weiss et ala 
and showed the commensurability oscillations in magne- 
toresistance for low magnetic fields and Shubnikov-de 
Haas (SdH) oscillations for high magnetic fields. Period- 
icity of the former ones can be well understood even in 
a semiclassical (SC) concept considering the drift of the 
cyclotron orbit centre in crossed electric and magnetic 
fields. This approach can also give gome quantitative 
predictions for the magnetorcsisiancea. An alternative 
formulation of the SC approacha relying on the break- 
down probability (tunnelling between two closed SC or- 
bits) oscillations is also possible. Gerhardts et ala di- 
agonalised the quantum-mechanical (QM) Hamiltonian 
(and employed the Kubo formula to compute conductiv- 
ity) finding the oscillating width of Landau bands (Lan- 
dau levels broadened by the weak modulation into nar- 
row cosine-like bands) to be the basic cause of the effect 
within the QM pictucp. QM approach was also applied 
by Vasilopoulos et alB. 

The SdH oscillations follow naturally from the QM 
concept owing to the quantization of free electron mo- 
tion in magnetic fields (Landau levels, LLs). This quan- 
tization has to be ad hoc assumed in the SC picture but 
once this is accepted, the SC theory provides a sufficient 
description in this case. 

Experiments on gated structures manufactured, by 
lithographic techniques performed by Bcton et alB al- 



lowed for investigating effects of stronger modulation. 
Compared to the previously mentioned experiments, the 
major effect was the quenching of the commensurabjlity 
oscillations. The SC approach^ was applicable agairfj, al- 
though quantum calculations^ were in better agreement 
with experiments. Later, the SC theories have beea ex- 
tended in order to account for anisotropic scattering^ and 
for modulation by magnetic fieldElEj (extending the older 
worksBO) . Ear QM approach to modulation by magnetic 
field see e.g.EE All together, the SC approach proved it- 
self to be suitable in these cases. 



More recently, samples with modulation potential due 
to MBE-grown structure were prepared (using cleaned 
edge overgrowth technique) by Deutschmann et alS3 In 
contrast to all former experiments, the modulation pe- 
riod d = 15 nm was shorter by almost one order of 
magnitude. Owing to this fact and also due to strong 
modulation (Fermi energy Ep ~ 4|i|, where t is the hop- 
ping integral between ground states in two neighbouring 
wells of the potential modulation), the lowest modula- 
tion miniband is well separated from higher minibands 
(condition for this is d < yj '3h 2 / \2\t\m)) . The magne- 
toresistance oscillations measured when the Fermi level 
lies between the modulation bands can be explained by 
no semiclassical model unless tunnelling betwecB. open 
trajectories is assumed (the breakdown formalisms men- 
tioned above is necessary at this place). It is thus ap- 
propriate to revert to a quantum mechanical description. 
Moreover, the miniband structure is simple now and al- 
lows thus for a good insight into the physics both on the 
SC and quantum-mechanical level. 
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bility of the SC approach will be shown to be hu> e g <C 2\t\ 
whereas w c ff = eB/m e g = eB j \jmm y . 

In the second part (Sec. |J) wc wm describe a fully- 
quantum— mechanical one-particle model and demon- 
strate that the gaps in the DOS emerging from this 
model coincide with extrema in the measured magnetore- 
sistance. Th en we will employ the linear response theory 
(section ID) in order to calculate the magnetoresistancc 
and we will compare it with the experiments. 



FIG. 1: 2DES with lateral periodical modulation and the 
effectively two-point-contact geometry. 

A. The System and Outline 

In this paper, we refer to experiments carried out 
on GaAlAs/GaAs structuropfo-et reported inl!3, see also 
more detailed description iiJia-ta. These are superlattices 
with strong unilateral short-period (electric-field) mod- 
ulation. The substantial difference to previous studies 
(see the Introduction) is the shortness of the modula- 
tion period: only the lowest modulation band (and not 
many of them) is occupied under those circumstances 
which makes the usage of quantum mechanics unevitablc. 
Moreover, absence of the higher modulation bands makes 
the model very transparent. We concentrate on magnc- 
toresistance measurements at different concentrations of 
electrons (which could be varied by a gate voltage over 
the range 0.5^5.0 x 10 11 cm -2 ). The system is sketched 
in Fig. |. 

The band structure in zero magnetic field calculated 
within the Kronig-Penney model can be well approxi- 
mated by 

E(k X ,ky) = ^~2\t\ COS kyd (1) 

with A\t\ « 3.8 meV and m equal to the effective mass 
of electrons in GaAs. The next modulation band is well 
above the Fermi level Ep for all accessible concentrations. 

If the Fermi level lies near the bottom of the band 
( — 2 1*| < E F < 2|t|), the system resembles a free 2DES 
(i.e. paraboloidal spectrum with modified effective mass 
in the modulation direction m y = h 2 /2\t\d 2 « 2.7m). If 
the Fermi level lies high above the modulation band edge 
(Ep 3> 2\t\) : the dispersion relation is like the one of an 
array of separated one-dimensional wires ( i.e. parabolic 
across k x and nearly constant along k y ) while the devi- 
ation from E = ti 2 k 2 ,/2m reflects the coupling of wires. 

The structure of this paper is the following. At the 
beginning we will review the semiclassical approach as a 
favourite tool used to describe magnetotransport exper- 
iments. We will also compare the zero-field density of 
states (DOS) with the B £ DOS computed in Sec. 0. 
It will give us illustrative examples of a situation when 
the SC theory is expected to be successful and of another 
situation when it should fail. The condition of applica- 



B. Semiclassical Approach 

The physical quantity of central importance in trans- 
port theories is the density of states g(E) (DOS) at the 
Fermi level Ep. It is known that its structure reflects 
features of the resistance (both as a function of B, for 
instance) but the relation between these two quantities 
is not simple. 

The SC theory attempts to explain the behaviour of 
electrons subject to a magnetic field in terms of the zero- 
field Fermi surfaces. We demonstrate that the zero-field 
DOS plus an extra quantization condition is a fairly good 
approximation to the realistic DOS (computed by our 
model, from Section |l|) at low magnetic fields (fiw e s <C 
2|t|). However, there is a drastic difference between these 
two densities of states for high magnetic fields indicating 
failure of the SC theory (see Fig. ||). 

Let us briefly review the SCI approach suggested by Lif- 
shitz and Onsager (see e.g.ll3). We construct the Fermi 
contour Ep = E(k x ,k y ) for a given Fermi level. The 
statements are that (1) the Fermi contour rotated by 
90 degrees and scaled by £ 2 = H/eB corresponds to the 
real-space trajectory of an electron (see Fig. ||) and (2) 
if the contour is closed, then it is allowed only if the 
magnetic flux passing through the area enclosed by the 
real-space trajectory is an integer multiple of the mag- 
netic flux quantum $o = h/e. In the QM picture, this 
quantization condition corresponds to the situation when 
Ep = hw c s(n + i) for some integer n, i.e. when the n th 
LL passes through the Fermi level. 
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FIG. 2: Real-space semiclassical trajectories of electron in 
magnetic field. - closed (Ef < 2|t|), 1 - critical (Ef = 2|t|), 
2 - open (Ef > 2|t|). The closed trajectories are elon- 
gated in the direction parallel to the wires by a constant fac- 
tor \/m y /m and due to non-parabolicity of the cosine band 
(Eq.II). 
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So as to be able to compare the SC and QM predic- 
tions let us now examine the densities of states. Note 
that by comparing the zero-field DOS of our system (for 
E -C 2|f|, i.e. in the region of quantized orbits and 2D- 
like behaviour) with the zero-field DOS of a free 2DES 
(g(E) = 2m e ft/Tth 2 including spin) we may deduce the 
modulation influenced effective mass m ff. This in turn 
determines the quantization condition Ep = huj e s(n + |) 
and thus all SC predictions can be made using the zero- 
field DOS of the system only. Based on the spectrum 
(Eq. pi) we can compute the zero-field DOS analytically 



9o(E) = { 



where 



4 / 2m 1 r- 
Pto"^ ,f>2|<1 



2m 



(2tt) 2 V h 2 \t\d 2 



(2) 



■K(y/l), -2\t\<E<2\t\ 



\ + E/2\t\ 



including the factor of two for spin (the dotted line in 
Fig. ||). K is the full elliptic function K(k) = F(§ , k) 



/; /2 (i 



; if^^dif, we recal]0 that K(0) = n/2, 



k sm* 
K(l) = oo. 

Let us focus on the weak-field case first (huj c s -C 2\t\) 
and discuss the influence of magnetic field B on the con- 
tinuous spectrum (Eq. p]), i.e. we try to estimate the 
DOS in magnetic fields without any calculation. On one 
hand, Landau levels (LLs) appear for Ep <C 2\t\ (2D- 
like paraboloid band structure, modified effective mass 
m y ). If the thermal energy is comparable to the Landau 
level separation (fc^T m tuv) the DOS becomes oscil- 
latory (leading to SdH oscillations) but approaches the 
zero-field DOS (see also Fig. ||). In other words, the 
oscillations missing in the zero-field DOS are exactly re- 
produced by the SC quantization condition (closed tra- 
jectories, Ep < 2\t\) which is fulfilled just when Ep lies 
in the middle between two LLs. 

On the other hand, narrow gaps open in the continu- 
ous spectrum for Ep » 2\t\ due to the slight corrugation 
(by the cosine term) of the almost ID-like parabolic- 
trough band structure. Numerical calculations (using 
model from Sec. Ill) show that these gaps are nar- 
row enough to disappear due to thermal broadening and 
the zero-field DOS matches the non-zero-field DOS per- 
fectly. The SC approach relying on non-quantized open 
trajectories (Ep > 2\t\) and predicting non-oscillatory 
magnetoresistance will therefore be successful once again. 

The picture is considerably different for strong mag- 
netic fields (huj e ff > 2\t\) at which the cyclotron radius 
R c approaches d. The numerically computed DOS shows 
no similarity to the zero-field DOS (Fig. ^ again). The 
non-zero-field DOS rather resembles g(E) cx 1/\Z~E of a 
single quantum wire (suppression of tunnelling between 
two neighbouring wires when classical cyclotron radii be- 
come comparable to d) with gaps both for E < 2\t\ and 
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FIG. 3: The density of states for B — (dotted) and ther- 
mally broadened (T w 1 K) density of states for low field 
(hu) e g <C 2\t\ or a <C 1, solid grey line) and high field (solid 
black line). 



E > 2\t\ owing to the tunnelling between the wires. 
Therefore we expect oscillatory magnetoresistance for 
both closed and open SC trajectories on contrary to the 
SC predictions (magnetic breakdown in the SC theory). 

The gaps occur at the boundary of the first magnetic 
Brillouin zone (1MBZ, see the section ||) and result from 
the periodicity of the dispersion relation in k x which 
in turn reflects the invaeiance of the QM Hamiltonian 
to magnetic translations^ (this is what remains from 
the full translational symmetry in the x direction after 
switching on the magnetic field). Magnetic breakdown 
can be included also in the SC picture through tunnelling 
between two open or two closed trajectoriesu (see Fig. |^). 
This is an ad hoc assumption though. 

The failure of the SC approach for strong magnetic 
fields is apparent also in another context. Once the Fermi 
level is set and the type of trajectory is determined, the 
same behaviour (either 2DES-like SdH oscillations or 
IDES-like no oscillations) is predicted for all magnetic 
fields. However, it is clear even on the SC level that if 
the cyclotron radius R c = kp£ 2 (kp is the Fermi wavevec- 
tor, E = ti 2 kp/2m) becomes comparable to d, tunnelling 
between the wires is suppressed and the system ought to 
switch to a quasi-lD regime. 

The QM approach is suitable for both low and strong 
magnetic fields. Moreover, it opens up the way to quan- 
titative calculations of the magnetoresistance and thus to 
data directly comparable to the experiments. 



II. MODEL 

Having chosen the Landau gauge A = (By, 0, 0) our 
system is described by the separable Hamiltonian (e = 

M) 



H=^-( Px + eByf + ^-p 2 y + V(y) , 
2m 2m " 



(3) 
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i.e. allowing to set ^(x, y) = ex.p(ikx)tp(y) for the eigen- 
functions. Our ansatz for the whole wavefunction is 



*(z,y) 



i 



exp(ikx) ^2 a 3 ( k ) <f(y ~ 3 d ) > ( 4 ) 



i.e. we use the ansatz \tp(k, n)} — J2j a j(^)\j) f° r V'(y) 
{n is the Landau index in the spectrum of (P) or (||) for a 
given k) where |j) is the ground state localized in the j— 
th well of the potential (corresponds to the Wannier state 
of the B = case). We have thus limited our model just 
to the lowest band in the modulation direction by this 
ansatz. The Fermi level lies always deep below the top 
of the modulation potential in our calculations. 

Next we use the tight-binding approximation (i.e. 
= = t5i_j±i, t < 0) and obtain the 

Hamiltonian matrix elements (see also Wulf et altB) 



Hjj — — 
3 2m 



K 2 ((k/K)+i) 2 5 lJ +tS 



K = d 



eB 



(5) 

In our model the real physical system is thus represented 
by the parameters t (hopping) and d (period) for the 
structure and of course B for the magnetic field. 

Note that (up to the scaling of k and energy) the prob- 
lem (^|) effectively depends only on the single parame- 
ter a: Eq. can be written in dimensionless form 



H tJ = \t\ a 2 ((k/K)+iyS, 



e 2 B 2 



2|t| 



\ 2\t\ 



i,j±l 
2 



(6) 



If we now assume the system to be infinite in the y direc- 
tion the matrix problem (j^) is mathematically equivalent 
to the Mathieu equation 

h 2 

- - ^0)2|t| cos Kx = Eibix) (7) 

2m 

i.e. ID Schrodingcr equation with potential W(x) = 
—2 |t| cos Xa— The detailed description of this model can 
be found inE2l. 

This allows us to determine the spectrum of (||) qual- 
itatively, and therefore also the DOS, even before we 
carry out the numerical calculations. Let the magnetic 
field be weak at first (a <C 1, see explanation below 
in this paragraph). Bands E(k) should appear (as a 
consequence of the Bloch theorem) and we may limit 
ourselves to 1MBZ, the first "magnetic Brillouin zone" 
k G {—hK, \K). For E <C 2|*| we expect those bands 
to be equidistantly spaced and narrow. The former fol- 
lows from W(x) w 2|*|(— 1 + \K 2 x 2 ) near the poten- 
tial minimum (the condition q«1 means just that the 
spacing of such states is <C 2|*|), the latter is due to the 
smallness of the overlap of low-lying states in two neigh- 
bouring wells of W(x). In contrast, for E 3> 2|*| we 
expect the spectrum to be almost like the one of ID free 
electrons E(k) = J2 i h 2 (k + iK) 2 /2m. The underlying 



non-constant potential V will manifest itself in the gaps 
which open at the Brillouin zone boundaries (k = zt^K). 

Let us now translate this analysis into terms of the 
original problem |^. Consider a fixed magnetic field (or 
constant a). For low energies we obtain nearly equidis- 
tant and sharp Landau levels (free 2D electron gas in 
magnetic field). At high energies we get a sum of almost 
ID densities of states, i.e. independent quantum wires 
(where B plays no role). The transition occurs around 
E = 2|f|. This agrees with the semiclassical theory ex- 
cept for the narrow gaps in the continuous quasi-free- 
lD-electron part of the spectrum. These indicate the 
magnetic breakdown for quantizing magnetic fields. 

Secondly, we focus on strong magnetic fields (a 1). 
On one hand, the states of Eq. ^ below 2|*| will become 
more widely spaced and even the lowest state will no 
longer have E <C 2|f|. We then expect even the lowest 
LB to be broad (solid black curve at Fig. ||) because its 
energy lies in an intermediate region between E 2|*| 
and E ^> 2|*|. On the other hand the states with E > 
2\t\ which approach the free electron states (and have a 
nearly free-electron parabolic spectrum) will have wider 
gaps (the smaller the lattice constant of a crystal, the 
wider the gaps). 

We may now compare the DOS being output from our 
model with the experimentally measured resistance, see 
Fig. |^. In the experiments, the Fermi energy was ad- 
justed by applying a gate voltage U g and we assume that 
Ep remains constant while magnetic field is swept. 

Experiments showed that the gate voltage is propor- 
tional to the concentration of electrons (N, see Sec. I_V). 
The proportionality constant is also approximatelly equal 
to the capacity of a parallel-plate capacitor correspond- 
ing to the gated structure. We took the concentration to 
be equal to the zero-field concentration and computed 
the corresponding Fermi energy, neglecting the localiza- 
tion effects. Looking at Fig. ^ we see the gaps in the DOS 
matching very well with the straight lines of magnctorc- 
sistance extrema which justifies this U g <-> Ep model of 
ours. 

We can see in the Fig. || that varying the magnetic 
field and keeping the Fermi energy (or U g ) constant and 
low the system behaves like an almost free 2D electron 
gas (sharp LLs, equidistant in 1/B) with modified ef- 
fective mass ^mniy (as predicted by SC theory). At 
intermediate Ep however, we observe that the narrow 
bands become broad at high B indicating an effective 2D 
to ID transition (when cyclotron radii become smaller 
than the modulation period). In contrast, the SC theory 
states that once the Fermi energy is set, a magnetic field 
cannot change the dimensionality of the system. When 
Ep is high the quantum wires are decoupled, but there 
still open gaps in the continuous spectrum (reflecting the 
tunneling between the open SC orbits). 

The lighter region around BkUT and U g f=a 0.45 V 
in the experimental data (Fig. ^ below) suggests that our 
model works no longer for strong magnetic fields (a 3> 
1) and low Fermi energies (even the lowest LB nearly 
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FIG. 4: Above: density of states (arbitrary units) in the 
tight-binding model for different magnetic fields and gate 
voltages. Below: measured resistance (logaritmic scale, ar- 
bitrary units). 



empty). This condition matches the situation when the 
cyclotron radius is much smaller than d and it means that 
our tight-binding model is inappropriate for magnetic 
fields strong enough to create LLs within one quantum 
wire (so that even the width of wires will be large enough 
for the electrons to behave as a 2DES inside one wire). 
These are, however, rather extreme conditions for the 
experiments shown in Fig. ^. 



III. TRANSPORT: KUBO FORMULA 

Knowing the eigenfunctions to the matrix problem |B| 
we pay use the linear response theory (Kubo formula, see 
e.g-Eil) to compute the conductivity tensor components 
(per system area A) 



o-u{E F ) 

Oxy{E) 



Tr 



irhe 2 
ihe 2 1 



Vi5 r (E F - H)vi5 r (E F - H) 



Tr 



v x G + (E)v v Sr{E-H) 



A 2 

v x Sr(E-H)vyG-{E) 



dN(E) 
' ~dB 



Here v x . y stand for the velocity operator components, 
5 r (E — H) = -(l/2wi)(G+(E) - G~(E)) and G ± (E) = 
(E — H ± iT). Isotropic scattering by a random impurity 
potential is taken into account by means of the complex- 
valued self-energy £ = A + iT and we neglect its real 
part. In the self-consistent Born approximation and as- 
suming T(E, B) to be small compared to the level sep- 
aration (typically fiw e ff) we can express components of 
<7 by means of the DOS g(E) and matrix elements of y 
(lu = eB/m) 



&xx(E) 



2 _ sgn(g(g)) 
irTd ' h ' g(E) 

.2 



AttT 2 e 2 ^ { (i/>(k,n')\y\iP(k,n)) 

— • (M • — • 9(E) ^ { E{ky) _ E ^ n) 



°yy(E) = ■ ^g(E) ]T (WM')ll#(M)> 



<Jxy(E) 



dN(E) 



dB 



(including spin degeneracy). The symbolic expression 
sgn(g(E)) I g(E) in the first term of a xx indicates that 
this term vanishes if g(E) = 0. Just to get simpler formu- 
lae, the additional assumption has been made that there 
are merely two points at the Fermi level E = E(k,n) 
within the,-J|MBZ. We reproduced the older results of 
Wulf et alS3 where the T-independent a xy (E) was cal- 
culated within the, same model but in a formally different 
way (Eqs. 5,6 incl). 

Note the structure of the components of a: in gen- 
eral, all the terms (except for the dN/dB term in 
u xy ) are products of the DOS (g) and some matrix 
elements of y. There are two contributions to the 
conduction parallel to the wires (a xx ): the first term 
(proportional to l/g) originates from the diagonal ma- 
trix elements S(k — k')hk/(eB) + (ip(k' ,n)\y\ip(k,n)) oc 
{ip(k',n)\v x \i;(k,n)) = (l/h)(dE/dk)S(k - k') cx l/g. It 
corresponds to the classical conductivity of a wire (open 
electron trajectories) or in other words it does not vanish 
owing to the unilateral modulation of the system: with- 
out modulation (free 2DES) the DOS comprises of delta- 
peaks (sharp LLs) and (sgng)/<? — > 0. For E < 2\t\, we 
obtain from our calculations LBs with non-zero width 
and this reflects tunnelling between two closed orbits in 
the SC picture (see Fig. ^), i.e. indicates a deviation 
from the quasi-2D behaviour. This contribution to con- 
ductivity is being suppressed by the impurity scattering 
(by 1/r) or in other words, persists even if there is no 
impurity scattering at all. 

The conduction perpendicular to the wires (a yy ) and 
the second term of a xx reflect the inter-LB transitions 
and appear thus due to impurity scattering (they are pro- 
portional to T) . These contributions might be viewed as 
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a consequence of tunnelling between the open SC orbits. 
There is no requirement that {^{n, fc)|y| , 0(m, k)) = for 
\n — m\ > 1 as in the limit of weak modulation, but we 
found these matrix elements to be decaying rapidly with 
growing \n — m\. In other words, inter-LB scattering 
occurs dominantly between neighbouring LBs. 

The Hall conductivity does not depend on the scatter- 
ing in the leading order at all. 



Magnetoresistance, high concentration 



Let us now cona 
tcring. Following). 18 '' 



depending only on energy (and B) 
consistent equation 



tc on the issue of impurity scat- 
!3 " we used an ansatz of self-energy 
This leads to a self- 



£(£) oc Tr 



1 



E — H — 



where the proportionality constant describes the strength 
of the impurity scattering. This equation yicldtoH the 
well- known result 



r 2 = —frio- 

2?r t 

for free 2DES (i.e. no modulation) and short-range scat- 
tering potential, r is the relaxation-time in the zero mag- 
netic field case (as in the Drude theory). In our calcu- 
lations we used this (i.e. T = jyB) as a phenomeno- 
logical ansatz which has already proven to be useful in 
explaining thc_magnetoresistance data by long-period 
superlatticesCjEH. Surprisingly enough, even this simple 
ansatz provides a very good qualitative agreement with 
the experimental data and makes the results of the cal- 
culations depend on the fitting parameter (7, scattering 
strength) in a very simple way. 

In order to obtain data comparable to experiments, we 
need to express the resistivity tensor components 



Byy 



Qxy 



(8) 




Magnetic field [T] 
Magnetoresistance, medium concentration 
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Magnetoresistance, low concentration 
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Magnetic field [T] 



' xx u yy 



■HI 



' xx w yy 



■nj 



A remarkable and important point of these formulae is 
that if scattering is weak (i.e. the inter-LB term oc T in 
a xx may be neglected), the denominators do not depend 
on r. 

Due to the particular sample geometry a two-point 
measurement is performed and thus the experimentally 
accessible quantity is some mixture of g yy and g xy . It 
seems plausible to assume that the voltage drop measured 
is a constant linear combination of the transversal and 
longitudinal voltage and thus the measured resistance is 

R CQyy ~\- Qxy 

with some dimensionless constant geometrical factor c. 
Thus if scattering is weak, the magnetoresistances to be 
compared to experimental data contain a single fitting 
parameter c/7. This comparison for different electron 
densities is shown at Fig. 0. 



FIG. 5: Magnetoresistance. Calculations and experiment. 
Electron concentrations (from the top) 5.1 x 10 11 cm -2 , 
3.8 x 10 11 cm" 2 , 2.65 x 10 11 cm" 2 . The lowest concentra- 
tion corresponds to Ef ~ 2\t\. 



IV. DISCUSSION 

Almost quantitative agreement is achieved for high 
concentrations (Fig. |5^.,b). Exact (no fitting parameter) 
match between the magnetoresistance extrema in theory 
and experiments (which has been emphasised already on 
the DOS level) shows that the tight-binding model cap- 
tures the essential physics in our experiments. 

Note also that the magnetoresistance approaches quan- 
tized values R = (\/v)h/e 2 for Ep lying in the gap. First, 
this allows to determine the concentration of electrons N 
in the system (compare filling factor of each plateau v 
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and its B). Second, it confirms our model of the relation 
between the gate voltage U g and Fermi level Ep. How- 
ever, as we may see best in Fig. ^p, there is still room 
for improving the model e.g. by taking the localization 
effects into account. 

The quadratic rise of resistivity at low magnetic fields 
due to the diagonal component (g yy ) has been pre- 
dicted within a SC pjpdclo and measured also at weak- 
modulation samplcsEH. Our results (with single fitting 
parameter at Fig. |5^,,b) deviate from this slightly be- 
cause of the form of r(£", B) oc \f~B. Assuming V inde- 
pendent on B in this region, the quadratic behaviour is 
reconciled. For intermediate fields the matrix elements 
of y start to play a role. 

Another feature of the tight-binding model is the ca- 
pability of explaining why the Hall plateaus show minima 
at high concentrations and maxima at low concentrations 
(Fig. ||c). Note, however, that even at this concentration 
E Ki 2\t\. Swapping the dominance of a xy and u xx yy 
(see esp. Eq. ||) is the cause for the qualitative change of 
the magnetorcsistance curve. Nevertheless, we are aware 
of the feebler match between theory and experiments at 
lower concentrations indicating that our ansatz for the 
self-energy is rather crude. 

V. CONCLUSION 

The cleaved-edge-overgrowth technology opens a new 
area of 2DES with short-period atomically precise uni- 
directional modulation for experimental studies. Such 



samples may be viewed as an array of coupled quan- 
tum wires which can be decoupled by applying reasonable 
(< 5 T) perpendicular magnetic field. From a theoretical 
point of view, these systems have a simple band structure 
and can be very well described by a simple tight-binding 
model. It was demonstrated that for weak magnetic field 
(hu> e g <C 2|i|) the SC approach argumenting with closed 
and open electron trajectories is expected to be fairly 
good and it is good indeed. On the other hand, this 
approach fails for quantizing magnetic fields (magnetic 
breakdown in the SC terminology). The fully quantum- 
mechanical tight-binding model offers both a good way 
to estimate the applicability of the SC approach and a 
reasonable description of the system for any magnetic 
field unless R c = kpl 2 <C d. This model is also very in- 
tuitive on the level of DOS analysis owing to the analogy 
with Mathieu equation (0). 

We found a good agreement between the experimental 
magnetorcsistance and calculations based on the linear 
response theory even with a very simple model for impu- 
rity scattering. We expect that the match between theory 
and experiment may be improved by internaly-consistent 
treatment of scattering due to randomly distributed im- 
purities. 
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